A new measurement of the properties of the rare decay K A 
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A large low-background sample of events (10300) has been collected for the rare decay of kaons 
in flight K + —* 7r + e + e~ by experiment E865 at the Brookhaven AGS. The decay products were 
accepted by a broad band high-resolution charged particle spectrometer with particle identifica- 
tion. The branching ratio (2.94 ± 0.05 (stat.) ± 0.13 (syst.) ± 0.05 (model)) x 10 -7 was determined 
normalizing to events from the decay chain K + —* tv + iy°; iy° — > e + e~'y. From the analysis of the 
decay distributions the vector nature of this decay is firmly established now, and limits on scalar 
and tensor contributions are deduced. From the e + e~ invariant mass distribution the decay form 
factor f(z) = /o(l + Sz) (z = M 2 e /rn 2 K ) is determined with S = 2.14 ± 0.13 ± 0.15. Chiral QCD 
perturbation theory predictions for the form factor are also tested, and terms beyond leading order 
0(p 4 ) are found to be important. 
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The decay K + — > 7r + e + e~ (K wee ) proceeds via a fla- 
vor changing neutral current and is highly suppressed 
by the GIM-mechanism. The decay rate was first cal- 
culated § assuming a short-distance s — ► d'j transition 
at the quark level. Later it was realized that the long- 
distance effects dominate the decay mechanism [||. Sev- 
eral recent calculations, which study the rate and invari- 
ant electron-positron mass (M ee ) distributions, were per- 
formed within the framework of chiral QCD perturbation 
theory (ChPT) , an approach which has been quite 
successful at modeling many decay modes of the light 
mesons 0. 

The first few events for this decay were observed at 
CERN §. Two subsequent experiments at the Brook- 
haven AGS, E777 § and E851 @, observed 500 and 
800 events, and measured branching ratios of 2.75 ± 0.26 
and 2.81 ±0.20 (xlO -7 ), respectively. We report here 
the results of a new measurement with larger acceptance 
and significantly increased statistics. Our data allow a 
detailed study of the decay form factor and a comparison 
with ChPT calculations and other models. 

The experiment was performed at the Brookhaven Na- 
tional Laboratory's AGS. The apparatus, a schematic 
drawing of which is shown in Fig. |l|, was constructed to 
search for the decay K+ — > 7r + ^ + e~ |lC^-|l3|]. It resided in 
a 6 GeV/c unseparated beam containing about 10 8 K + 
and 2 x 10 9 7r + and protons per 1.6 second AGS pulse. 
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FIG. 1. Plan view of the E865 detector. PI - P4: propor- 
tional chambers; CI, C2: Cerenkov counters; ABCD: scin- 
tillator hodoscopes. The beam passes through holes in the 
calorimeter and muon stack, dead regions in the proportional 
chambers, He bags and H2 filled beam tubes in the Cerenkov 
counters, which are not shown here. 

Downstream of a 5m-long evacuated decay volume a 
dipole magnet separated the trajectories by charge, with 
negative particles going mainly to the left. This was fol- 
lowed by a spectrometer consisting of four proportional 
chambers surrounding a dipole magnet with a 0.833 Tm 
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field integral, which determined the momenta and tra- 
jectories of the decay products. Particle identification 
was accomplished with Cerenkov counters filled with H2 
on the left (C1L and C2L), and CH 4 on the right (C1R 
and C2R) at atmospheric pressure; with an electromag- 
netic calorimeter of the Shashlik design Q, consisting 
of of 600 modules, each 11.4 cm by 11.4 cm by 15 ra- 
diation lengths in depth, arrayed 30 horizontally and 20 
vertically; and with a muon range stack consisting of 24 
planes of proportional tubes situated between iron plates. 
With these components, electrons were identified as hav- 
ing light in the appropriate Cerenkov counters and energy 
in the calorimeter consistent with the measured momen- 
tum of the trajectory. Pions were identified as having 
no light in the Cerenkov counters, and an energy loss in 
the calorimeter consistent with that of a minimum ioniz- 
ing particle or a hadron shower. Compared to the most 
recent experiment the E865 spectrometer had im- 
proved particle identification capabilities, spectrometer 
resolution and larger, more uniform acceptance. 
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FIG. 2. M ee distribution for accepted Kdai events with 
(lower curve) and without (upper curve) high-mass trigger 
required. The Monte Carlo simulation (histogram) includes 
events from other K + decays with a ir° in the final state. 

The first-level trigger for the experiment is based on 
three charged particle hits using hodoscopes D and A 
(arrays of 10 and 15 slats each on either side of the beam 
line, respectively), and the calorimeter. Most of the rate 
at this stage comes from accidentals, and K + — > 7r + 7r + 7r~ 
(K T ) decays. For K„ ee decays the next level trigger re- 
quired Cerenkov counter signals on each side of the detec- 
tor. This trigger was dominated by events from the decay 
chain K + — » 7r + 7r°; it — > e + e~7 (Kdai) with low invari- 
ant mass M ee . To enhance the K vee fraction a high mass 
trigger was configured, in which events with small ver- 
tical separation of calorimeter hits were prescaled. This 
reduced the number of Kdai triggers, while keeping 85 % 



of the high M ee events. Figure |^ shows the invariant 
mass spectra for the prescaled monitor Kdai events and 
those, which were accepted by this trigger and were used 
for normalisation. 
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FIG. 3. Scatter plot M ee versus M^ ee for Knee candidates. 
Insert: M neB mass for K^ee candidates with M ee > 0.15 GeV. 
The histogram shows the Monte Carlo simulation. 

In the analysis, the basis for the selection of both nor- 
malisation and signal events was the unambiguous iden- 
tification of a positive pion and lepton pair with tra- 
jectories from a common vertex located within the de- 
cay volume. With a cut M ee > 0.15 GeV we have a 
nearly pure sample of K vee events with a three particle 
invariant mass M wee « rriK, as shown in Fig. [3J. The 
onset of the Kdai events can be seen for M ee < 0.13 
GeV and M^ ee < mx, the latter because of a missing 
photon. For K wee events, the reconstructed K + was re- 
quired to come from the production target within the 
limits inferred from studying K T decays. Our final sig- 
nal sample contains 10300 K^ee candidates including 1.2 
% background events. Our normalisation sample, after 
prescaling, contains 10 5 Kdai candidates including a 17 % 
contribution from two other K + decays giving 7r°-Dalitz 
pairs (K + — ► 7r°/i + i^ with the [i + treated as a 7r + , and 
K+ -> 7T + 7r 7r°). The insert in Fig. | exhibits a M nee 
mass resolution of a = 5.7 MeV, in good agreement with 
the Monte Carlo simulation. The calculated M ee resolu- 
tion is a — 4.8 MeV nearly independent of M ee . 

The acceptance for K vee and Kdai events was deter- 
mined to be 0.73 % and 0.85 %, respectively, with a 
Monte Carlo simulation, which included the geometry 
of beam line and spectrometer, and the separately mea- 
sured efficiencies and responses of scintillators, propor- 
tional chambers, Cerenkov and shower counters. The in- 
fluence of the efficiencies on the acceptance nearly cancels 
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in the normalisation, since K nee and K^ai events contain 
the same final states, and their different spatial distri- 
butions produce only minor differences, and are taken 
into account. Figures || and |^ are examples of the many 
control plots which attest the excellent quality of the sim- 
ulation. For Kdai events the matrix elements of |l6| were 
used for ir° — > e + e~7, and the theoretical input to the 
Knee simulation is discussed below. Figure |2| also demon- 
strates that the major difference between K ve e and Kdai 
events, the high-mass trigger, is correctly accounted for, 
in magnitude as well as in shape. 
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FIG. 4. Angular (left) and invariant mass (right) distribu- 
tions for K^ee events (data points) compared to Monte Carlo 
simulations (histogram) assuming a pure vector interaction. 
A linear form factor parametrisation with S = 2.14 is used. 
The dashed histogram (right) corresponds to a constant form 
factor (5 = 0). 

The essential distributions necessary for the interpre- 
tation of our data are shown in Fig. [|. Since the decay is 
supposed to proceed through one photon exchange, i.e. 
by a vector interaction (V) with a decay amplitude || 
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expects an angular distribution proportional to 
sin" 9, where 9 is the angle between the positron and 
pion momentum vectors in the center of mass of the e + e~ 
pair. The presence of other decay mechanisms, however, 
may produce small admixtures of scalar (S) or tensor (T) 
terms 17 . The corresponding decay amplitudes 
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lead to either a constant (S) angular distribution or one 
proportional to cos 2 9 (T). Here Gp is the Fermi constant, 
P = PK+Pin 1 = Pk — ?V> and the form factors fv.s,T are 
dimensionless functions of z — M 2 e /m 2 K . Fitting a two- 
dimensional distribution whose projections are shown in 



Fig. ||, we find good agreement with a vector interaction. 
At 90 % confidence level at most 2 % of the branching 
ratio could result from either scalar or tensor interaction 
corresponding to constraints \fs\ < 6.6 • 10~ 5 and |/t| < 
3.7- 1(T 4 . 

Neglecting terms proportional to m 2 /M 2 e the mass dis- 
tribution for a vector interaction can be described by 
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where i](a, b, c) = a 2 + b 2 + c 2 — 2ab~2ac — 2bc. The form 
factor fv{z) can be determined by fitting the observed 
spectrum in the experimentally accessible range 0.1 < 
z < 0.51. We have used two different parametrisations 
of the form factor, one model independent (Eq. ||) and 
the other derived from ChPT §] (Eq. |): 

8'z 2 ) , 



f v (z)=Ml + Sz 
fv(z) = a+ + b+z + w™(z) 

f ,5, 5',a + ,b 



(2) 
(3) 

and w 7 ' 11 is the 



are free parameters |18 
contribution from a pion loop graph given in 

Figure |E| displays the form factor, which is extracted 
from the ratio of mass distributions for the measured 
events to events simulated with a constant form factor. 
The results of our fit to the linear and ChPT ansatz are 
superimposed. Figure |] shows the spectrum itself. The 
parameters and branching ratios are given in Table §. 
The different contributions to the systematic uncertainty 
of our results are listed in in Table H. Radiative correc- 
tions have been included in the simulation following [ p"5| . 
This increased the branching ratio by 5.5 % and the lin- 
ear slope by 4 %. 
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Eq. 



X 2 /n.do 



fo 0.533 ±0.012 
S 2.14 ±0.13 
6' 

22.9/18 



0.591 ±0.027 a+ 
0.97 ±0.44 b+ 
1.99 ±0.67 
16.6/17 



-0.587 ±0.010 
-0.655 ± 0.044 

13.3/18 



BR (total) 2.884 ± 0.037 2.991 ± 0.058 2.988 ± 0.040 
BR (meas.) 2.015 ± 0.020 

TABLE I. Summary of results for the fits to the measured 
form factor. The branching ratios (BR) are given in units of 
(10~ 7 ). BR (meas.) corresponds to M ee > 0.15 GeV and is 
within the statistical error independent of the parametrisation 
of the form factor. 



Source 


obrIBR (%) 


a s /5 (%) 


Normalisation to ty u — > e + e~7 


2.8 




Radiative corrections 


1.4 


0.5 


Reconstruction 


3.0 


6.5 


Background subtraction 


1.5 


2.3 


Total 


4.5 


6.9 



TABLE II. Summary of systematic errors. The reconstruc- 
tion and background subtraction error result from adding sev- 
eral smaller contributions linearly. The total error is the sum 
of the four groups in quadrature. 
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The inclusion of the quadratic term in the model inde- 
pendent ansatz improves the quality of the fit. The values 
of S and 6' are however strongly correlated. Our result 
5 = 2.14±0.13±0.15is in fair agreement with the less pre- 
cise results of experiment E777 6 = 1.31±0.44. The 
large value of 5 is in contradiction with the meson dom- 
inance class of models in which the form factor 
slope is given in first approximation as S = (mx /m p ) 2 . 

The ChPT parametrisation of Eq. || includes the con- 
tribution from the pion loop graph. This term has a 
strong z dependence, which has not been calculated ex- 
plicitly in most earlier models. In fact, in the first, lowest 
non-trivial order ChPT calculation (0(p 4 ) 0) this term 
was assumed to be the only significant source of z depen- 
dence of the form factor (i.e. &_|_ = 0). Our data show 
that this is a poor approximation. Significant z depen- 
dence from other terms is expected in next-to-leading- 
order ChPT §||. The amplitude linear in z (Eq. |) 
is thought to represent all contributions other than the 
pion loop term, which ChPT cannot calculate ||. The 
substantial reduction in the value of % 2 which is observed 
when the pion loop term is included provides direct ex- 
perimental evidence for its small but important contri- 
bution. 



1.4 



~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



> 



1.2 - 



0.8 



0.6 



0.4 




0.1 



0.2 



0.3 



0.4 



Z = (M ee /M K ; 



0.5 

,2 



FIG. 5. The measured form factor |/v(«)| versus z. The 
solid line shows the linear fit (Eq. |^ with 8' = 0), the dashed 
line the next-to- leading-order ChPT fit (Eq. [|. 

Summarizing the results of our analysis we conclude: 
(i) the experimental data are consistent with a vector 
model for the interaction, (ii) the slope of the form factor 
is significantly larger than meson dominance or leading 
order ChPT models predict, and (iii) although a linear 
approximation of the form factor is reasonable, our data 
indicate a nonlinearity of the form factor which is fit well 
by the ChPT loop term. Our final result for the to- 



tal branching ratio is (2.94 ± 0.05 (stat.) ± 0.13 (syst.) ± 
0.05 (model)) x 10~ 7 . Here we include the model depen- 
dence of the extrapolation into the low mass region not 
covered by our detector by taking the average of the ex- 
treme values in Table |[ 
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